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ارزیابی تبخیرتعرق مرجع با استفاده از روش‌های داده‌کاوی و مقایسه آن با نتایچ سامانه 
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آزاده صداقت *- نیازعلی ابراهیمی پاک "- آرش تافته ۲ سیده نر گس حسینی* 
تاریخ دریافت: ۱۴۰۱/۰۶/۱۴ 


تاریخ پذیرش: ۱۴۰۱/۰۹/۱۵ 


چکیده 


برآورد دقیق تبخیرتعرق مرجع (۲0) یکی از عوامل مهم برای محاسبة نیاز آبی و آب مصرفی گیاهان زراعی و باغی است. پیچیدگی فرآیند 
کر ضرق و زاگ ی به‌دامم‌های ,سای برازه قق این مقر راشای وه اس فیرگی خی عم قطیی داتی و یار به انااشت ات 


قلیمی متنوع در برآورد 0 از دلایلی بوده‌اند که باعث شده پژوهشگران به‌سوی روش‌های داده‌کاوی همچون شبکه عصبی مصنوعی (۸((5) 
جنگل تصادفی (1۳) و ماشین بردار پشتیبان (8۷[۷1) روی آورند. در این تحقیق, داده‌های هواشناسی در بازه زمانی ده ساله (۱۳۸-۱۳۹۹) از 
پستگاه‌های هواشناسی استان قزوین جمع‌آوری شد. ابتدا مقادیر 31 در سامانه نیاز آب که از روش پنمن-مانتیث محاسبه شد. استخراج گردید. سپس 
ین مقادیر به‌عنوان مقادیر واقعی (اندازه‌گیری شده) با مقادیر تخمینی بدست آمده با روش‌های داده‌کاوی (۸(5 3۳ و )٩۷1۷‏ ارزیابی شد. جهت 
عتبارسنجی نتایج بدست آمده» داده‌های هر ایستگاه به دو مجموعه آموزش (دوسوم داده‌ها) و آزمون (یک‌سوم داده‌ها) تقسیم شدند. نتایج بررسی‌های 
آماری و دياگرام نشان دادنه در هر سه روش استفاده شده با در نظر گرفتن تمامی پارامترهای هواشناسی (میانگین دمای هواء ميانگین رطوبت نسبی, 


ساعت آفتابی و سرعت باد) به‌عنوان ورودی مدل. در ایستگاه سینوپتیک قزوین و ایستگاه کلیماتولوژی نیروگاه رجایی, در هر دو مرحله آموزش و آزمون» 
60 با دقت بالاتری برآورد شد. همچنین در این تحقیق دقت نتایج روش ۸5 نسبت به دو روش دیگر به‌طور نسبی بالاتر بوده است. در هر دو 
مرحله آموزش و آزمون مقادیر ۷۲ و 12 بداست آمده از روش 5( در ایستگاه سینوپتیک قزوین برابر و به ترتیب برابر ۰/۲۱ و ۷« و در 
ایستگاه کلیماتولوژی نیروگاه رجایی برابر و به‌ترتیب برابر ۰/۱۰ و ۰/۹۷ می‌باشد. به‌طورکلی نتایج نشان داد که میانگین دمای هوای روزانه مهمترین 


پارامتر هواشناسی تأثیرگذار در برآورد 7 می‌باشد. 
واژه‌های کلیدی: تبخیرتعرق مرجع. داده‌کاوی سامانه نیاز آب 
مقدمه 


تبخیرتعرق فرآیند از دست رفتن آب از خاک و سطوح محصول به 
اتمسفر توسط فرایندهای تبخیر و تعرق است. تبخیرتعرق یکی از 


۸ ۲ ۲ و ۴- به‌ترتیب محقق, دانشیا, استادیار و محقق بخش مدیریت آب در 

مزرعه موسسه تحقیقات خاک و آب سازمان تحقیقات» آموزش و ترویج کشاورزی 

کرج» ایران 

(#- نویسنده مسئول: 0211.6017ع ۵ ۸2۵0۳86022102665 :۵21ظ۴ظ) 
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مهم‌ترین اجزای چرخه آب است و عامل کلیدی برای کشاورزی» 
برنامه‌ریزی آبیاری و منابع آب است. پیش‌بینی دقیق و سریع 110به 
تجزیه وتحلیل تغییرات محیطی کمک می کند و برای ابیاری گیاه 
ارزیابی محیطزیست در مقیاس‌های مختلف و هیدرولوژیکی و مدل 
اکوسیستمی ضروری و بسیار مهم است (2021 ,طع1؟" 240 عظع۳). 
پیچیدگی فرآیند تبخیرتعرق و وابستگی آن به متفیرهای هواشناسی» 
عدم دسترسی به تمامی داده‌های هواشناسی و عدم تعمیم‌پذیری یک 


۴ نشریه آب و خاک. جلد ۰۳۶ شماره ۶ بهمن - اسفند ۱۴۰۱ 


است. به‌طورکلی برآورد 770 را می‌توان از طریق روش‌های مستقیم و 
غيرمستقيم انجام داد (2019 ,.۵1 2۶ ع118). ۳0 را می‌توان به‌طور 
مستقیم بر اساس روش‌های تراز انرژی و انتقال جرم بخار آب 
اندازه‌گیری کرد بااین حال, هزینه زیادی دارد. روش‌های غیرمستقيم 
برای تخمین ۳0 با توسعه معادلات تجربی مانند-۵۳۵0ع۳ 
(۳(۷) طانعا10 آغاز شد (1965 ,طاه۷]0۸1). بر اساس روش 
۲ ۲۸056 یک روش استاندارد برای تخمین 310 منتشر کرد. 
روش ۳۲ ارائهشده در 1۸056 توسط بسیاری از محققین به‌عنوان 
روش پایه در محاسبه ۳0 استفاده شده است معن12 20 (عع2). 
(2019 ,۸ ۰۶ تتنطگ ,2019 روش (۳ علاوه بر دقت بالا در 
تخمین 10 دارای محدودیت‌هایی مانند نیاز داشتن به پارامترهای 
آب و هوایی زیاد می‌باشد. از طرفی ایستگاه هواشناسی کمی وجود 
دارد که همه این پارامترها را باهم ارائه دهد ,۵ ۶ ۲اعوعانقد۳) 
(2014. برای غلبه بر این مشکل, بسیاری از محققان مدل‌های 
تجربی را توسعه داده‌اند که نیازهای کمتری به داده‌های آب و هوایی 
دارد. 
رویکردهای دیگری که اخیراً توجه محققین را به خود جلب کرده 
است» روش‌های مبنی بر هوش مصنوعی مانند شبکه عصبی 
مصنوعی (۱)۸(۷(۲5 (2018 ,.۵1 6۶ 02۷31)» ماشین بردار پشتیبان 
۲۵6۷3 (2013 .1 ۰۶ ت12 و درخت‌های تصمیم‌گیری مانند 
(۲)۴۳ (2020 .۰1 21 132272) می‌باشد که کاربردهای زیادی در 
تخمین 1710داشته‌ان. با طلهور فناوری یادگیری ماشین و هوش 
مصنوعیء محققان شروع به کشف چگونگی ترکیب الگوریتم های 
هوشمند با روش های تخمین سنتی برای تخمین دقیق تر و موثرتر 
0 کردند. یونس و همکاران (2020 .۵1 ۶ وعطنا» 710 روزانه را 
بر اساس شرایط آب و هوایی با استفاده از معادلات تجربی و 
روش‌های مختلف داده‌کاوی پیش‌بینی کردند. با توجه به نتایج آن‌ها؛ 
معادله مبتنی برتابش بهتر از سایر معادلات تجربی و بالاترین ضریب 
همبستگی برای روش شبکه عصبی و حداقل خطاها برای تابع مبنای 
شعاعی ٩۷1۷۲‏ به‌دست‌آمده است. چن و همکاران ..۵1 6 عع) 
(2020 ۲0 روزانه را بر اساس داده‌های هواشناسی محدود با 
استفاده از سه روش بادگیری عمیق» دو روش یادگیری ماشین و هفت 
معادله تجربی برآورد کردند. نتایج آن‌ها نشان داد که وقتی 
ویژگی های مبتنی بر دما در دسترس باشد» مدل های یادگیری عمیسق 
به‌طور قابل توجهی بهتر از مدل‌های تجربی مبتنی بر دما عمل می- 
کنند. دوس سانتوس فاریاس و همکاران ( 6 عدتتد۲ دماصدگ دم 
0 ,.21) عملکرد روش‌های یادگیری ماشین برای تخمین 1710 
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روزانه در کشور برزیل ارزیابی کردند. نتایج آن‌ها نشان داد که 
روش‌های یادگیری ماشین» حتی در غیاب برخی متغیرهاء 110را با 
دقت بالا برآورد می‌کنند. متار (2018 ,۷]21627) هشت مدل برای 
تخمین 170 با استفاده از الگوریتم ژنتیک ایجاد کرد. نتایج این 
تحقیق نشان داد نتایج برآورد بسیار نزدیک به تخمین با مدل1۸0- 
۳[۷]۲-6 بوده است. 

استفاده از روش های هوش مصنوعی در تخمین ۳0 با دقت بالا 
در سال‌های اخیر رواج یافته است اما پیچیدگی این روش‌ها و نیاز 
آن‌ها به داده‌های هواشناسی بیشتر کاربرد آن‌ها را برای مناطق 
مختلف دشوار می‌کند (2021 م12 200 ۳۲۵۸۵). بنابراین» در مقالة 
حاضر 1 از چهار متغیر هواشناسی (میانگین دمای هوای («:1)* 
ميانگین رطوبت نسبی (میب[۵8 سرعت وزش باد در ارتفاع دو 
متری (۷۷۰5) * و ساعت آفتابی (:/:۲)5) که از ده ایستگاه هواشناسی 
سینوپتیک و پنج ایستگاه کلیماتولوژی استان قزوین به‌دست‌آمده 
است» با استفاده از روش‌های داده‌کاوی تخمین زده شد. هدف از این 
مطالعه نشان دادن شایستگی روش‌های داده‌کاوی جهت رسیدن به 
یک مدل جامع و ساده, و معتبر با حداقل داده‌های آب و هوایی برای 
تخمین ۲0 روزانه می‌باشد. 


مواد و روش‌ها 
موقعیت محدوده موردمطالعه 

استان قزوین در حوزه مرکزی ایران بین ۴۸ درجه و ۴۵ دقيقه تا 
۰ درجه و ۵۱ دقیقه طول شرقی از نصف‌النهار گرینویچ و ۳۵ درجه و 
۴ دقیقه تا ۲۶ درجه و ۳۸ دقیقه عرض شمالی نسبت به خط استوا 
قرار دارد (شکل ۱). استان قزوین یکی از استان‌های مرتفع ایران است 
و ارتفاع بیش از ٩۰‏ درصد سطح استان بیش از متوسط ارتفاع 
(۱۲۰۰متر) است. اختلاف شدید ارتفاع و وجود ناهمواری از عوامل 
اصلی تنوع آب‌وهوای استان قزوین است. بیش از نیمی از استان 
قزوین یعنی ۵۲ درصد دارای اقلیم خشک فراسرد تا سرد و ۱۸ درصد 
از سطح استان دارای اقلیم مدیترانه‌ای فراسرد تا سرد می‌باشد. در این 
پژوهش از داده‌های مربوط به دوره ده سالةّ (۱۳۹۹-۱۳۸۹) استفاده 
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شکل ۱- موقعیت جغرافیایی استان قزوین 
۸۵ 027۷11 0۴ ممناه)م۱ 066۵۵۵ -1 ۲12۱۲۵ 


برآورد مقداره[1 سامانه نیاز آب گیاهان 

مقدار 1710 در این پژوهش از بانک اطلاعاتی سامانه نیاز آب 
موسسه تحقیقات خاک و آب کشور که به‌منظور تعیین نیاز آبی 
استاندارد و واقعی گیاهان باغی و زراعی طراحی‌شده است. استفاده شد 
(۸)۵://۳:۷۷۲۰1۲ظ). استفاده از این سامانه به‌این‌علت می‌باشد که در 
حال حاضر دانش فنی جامعی که بتواند از بانک‌های اطلاعات کلان 
کشور برای برآورد و تعیین نیاز آبی» آب مصرفی و برنامه‌ریزی آبیاری 
در مناطق مختلف» استفاده شود» وجود ندارد. این سامانه بانک 
اطلاعاتی خوبی می‌باشد که در زمان کم» حجم زیادی از اطلاعات 
منطقه را در اختیار کاربران قرار می‌دهد. مقدار ۳0 استخراج‌شده از 
سامانه نیاز آب بر اساس روش پنمنمانتیث فائو ۱۵۶و در مقیاس 
زمانی روزانه محاسبه‌شده است» که به‌عنوان مقدار واقعی (اندازه‌گیری 
شده) در مدل‌سازی با روش‌های داده‌کاوی در نظر گرفته شد. 


توصیفی از مراحل تخمین 1710و روش‌های مدل‌سازی 
برای ایجاد مدل در تمام ایستگاه‌های هواشناسی» هم‌راستایی بین 
متفیرهای مستقل (فرض صحت رگرسیون خطی) از طریق آماره 
عامل تورم واریانس (۲)۷۳ و با استفاده از نرم‌افزار با 8۴55.26 
بررسی شد. متغیرهایی که ۷1۳ آن‌ها کمتر از ۱۰ باشد مشکل 
هم‌راستایی نخواهند داشت (2013 ,ع۳10010). جهت اعتبارسنجی 


1- ۲۸۵-56 ۳۵۴۵۵۴ 
2- ۱۷ 2۳121066 ۲۴120100 ۲ 


نتایج به دست‌آمده از این تحقیق. داده‌های هر ایستگاه به دو مجموعه 
آموزش (دوسوم داده‌ها) و آزمون (یک‌سوم داده‌ها) تقسیم شدند. 

در این پژوهش از روش‌های شبکه عصبی مصنوعی (۸((5)» 
ماشین بردار پشتیبان (5۷۸) و جنگل تصادفی (2۳) جهست 
مدل‌سازی ۳۲0 استفاده شد. در این پژوهش برای مدل‌سازی ۳10 با 
روش‌های ۸۱۷۶ 5۷/۲ و 13۳ از نرم‌اف زار 5۲۸۲16۲16۸12 
استفاده شد. شبکه عصبی پرسپترون چند لایه (۷11۴)" و تابع پایه‌ی 
شعاعی (3۳)" دو نوع شبکه از روش ۸۲ می‌باشد که در این 
پژوهش مورد استفاده قرار گرفت ( ,طماعزط :1996 ,.۵1 ۶ 20ع۳۱۸ 
0 روش ۸۸ در ۱۰۰ تکرار جت کاهش خطادر 
مدل‌سازی» تنظیم شد. میانگین ۵۰ خروجی از ۱۰۰ تکرار که خطای 
کمتری دارنده به‌عنوان خروجی نهایی برای تخمین 1710 انتخاب شد. 
روش ٩۷۷۲‏ دارای عامل‌های تنظيم‌کننده شامل ) و 1 هسته 
می‌باشد که دارای چهار نوع کرنل شامل کرنل خطی, چندجمله‌ای, 
شعاعی و حلقوی است (2000 ,۵0:11 ۷). روش 18۳ یک روش جدید 
و قدرتمند است که پیشرفت‌های فراوانی را در فنآوری داده‌کاوی 
داشت. اساس این روش توسط (2001) ,ص2حصع:ظ ارائه شده است. 


معبارهای ارزیایی دقت و قابلیت اعتماد مدل‌ها 
برای بررسی دقت و قابلیت اعتماد مدل‌های ایجادشده از 
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معیارهای آماره‌های ضریب تعیین (*1۳)" مجذور میانگین مربعات خطا 


("115۳» و میانگین انحراف خطای (0۷19۳۳» مجذور میانگین 
خطای نرمال شده (08۷18۳) و ضریب کارایی مدل (۳۳) 
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از لحاظ آماری حداقل مقدار 120۷15 زمانی مقادیر تخمینی با 
مقادیر اندازه‌گیری برابر صفر باشد مدل دارای دقت بیشتری خواهد 
بود (2020 ,.41 61 ۰)۸20 ۷1۳913 عمدتاً برای نان دادن اریب در 
مدل, استفاده می‌شود. وقتی مقدار آن مثبت باشد نشان‌دهندة بیش 
برآوردی و مقدار منفی آن نشان‌دهندة کم‌برآوردی تخمین می‌باشد 
(1988 م,حعطنععنومن). از ۱۳۴/8۳ به‌عنوان شاخص خطای 
اظر ی ال اسفاهمتس شیاف اک نار ان کت ال ۱۶ ۸ افتهه 
نتیجه تخمین عالیء بین ۱۰ تا ۰/۲۰ خوب؛ بین ۲۰ تا ۰/۳۰ قابل 
قبول؛ و برای مقادیر بیشتر از //۲۰ ضعیف می‌باشد ( 4مد «طعح5 
0 ,1۳۳21). شاخص » شاخص توافق ویلموت می‌باشد که از 
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دامنه صفر (بدون توافق) تا یک (توافق کامل) می‌باشد که میزان 
درستی و صحت برآوردها را نشان می‌دهد ( ,۵ 6 نامصال ۱۷ 
2 همچنین مقدار1۳ از منفی بی‌نهایت تا یک به ترتیب برای 
بدترین و بهترین برآورد تغییرمی‌کند و کفایت و کارآیی مدل را نشان 
می‌دهد و بالای ۰/۵ عملکرد قابل‌قبول مدل را نشان می‌دهد 
(2020 ,تماعلماتصدظ 4صه عامحتدنتمک). 


توصیف آماری 

ویژگی‌های آماری متفیرهای موردمطالعه در همه ایستگاه‌های 
سینوپتیک و کلیماتولوژی در جدول ۱ آورده شده است. تمامی 
داده‌های ورودی و خروجی در همه ایستگاه‌ها در این تحقیق. دارای 
دامنه تغیبرات بالایی (اختلاف بین کمترین و بیشترین مقدار) می- 
باشند» که نشان‌دهندة این است که باعث افزایش قابلیت تعمیم‌پذیری 
مدل می‌شود (جدول 4۱ به‌منظور نشان دادن پراکندگی نسبی داده‌ها 
به میانگین» از ضریب تغییرات (۷:))استفاده شد. همبستگی بین 
عامل‌های ورودی مدل با ۳10 بررسی شد (جدول ۲). نتایج 
همبستگی نشان داد. 710 در همه ایستگاه‌ها با ساعت آفتابی» سرعت 
باه ميانگین دمای هوا همیستگی ثبت و معنی‌دا, و با میانگین 
رطوبت نسبی همبستگی منفی و معنی‌داری دارد (جدول ۲) با افزایش 
تابش خورشید و بالا رفتن دمای هوا و هم‌چنین افزایش سرعت باد 
مقدار ۳10 افزایش‌یافته و با افزايش رطوبت هوا مقدار ۳10 کاهش 
می‌یابد (2020 ,.۵1 6۶ 6هذل). 


ویژگی‌های روش مورداستفاده در این پژوهش 

شبکه عصبی پرسپترون چند لایه (۷11۳) و تابع پایه‌ی شعاعی 
(62۲» دو نوع شبکه از روش ۸۷5 می‌باشد که در این پژوهش 
در ایجاد مدل مورداستفاده قرار گرفت (2000 ,۳16608). در نتایج 
بدست آمده از روش ۸۵ در همةایستگاه‌هاء تعداد شبکه 1۳ 
استفاده‌شده بیشتر از شبکه 8۳ می‌باشد. بنابراین شبکه ۸11.۳ 
نسبت به 188 در برآورد مقدار 1710 از عملکرد بهتری برخوردار 
است. در روش ۸5 از توابع فعال‌سازی مختلفی در لایه‌های میانی 
و خروجی استفاده شد. که ۲0۳6/11 و 1200 تابع فعال‌سازی 
غالب در لایه پنهان و 10000117 تابع فعال‌سازی غالب در لایه 
خروجی بوده و بهترین تابع برای توصیف روابط بین متغیرهای ورودی 
و 1770 می‌باشند (جدول ۲). یکی از مهم‌ترین ویژزگی روش 16۳ 
ساده‌سازی مدل از طریق تعیین اهمیت عامل " (تخمین‌گرها) می‌باشد 
(جدول ۴). 
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ساعت ."سرت ..."لگ ميانگیر ساعت._. سرعت . کت ميانگید 
۳ باد رطوبت ِِ ایستگاه‌ها ۳ باد رطوبت ث_ ایستگاه‌ها 
آفتابن ۱ ۱ و هو عصمناداه افتابی ۱ 1 دمای هو عصمناعاه 
اس ۲۷5 )سس «اسه. ۲۷5 0 
(1 _ (66) «ب«1۳ (1 _ (66) «ب«1 
۲ 1 
وی و0 060 09 کوهین ۱0 69 0- و09 بیک ۳10 
صرصرصرصپصپصپصپصپصپصپپپ ‏ 0 پپبب ۹ | 
اسماعیل آباد 
ک ا ‏ #ک04 05 09 معلم کلایه م04 060 و0 ممموط  .‏ 80 
ههام صصعل۷]۵2[ ۳-۷ 
۱0 02 و09 لبرز 0 011۴ 072 *#وع0 آوج ۳10 
سس 0017 رب ۱9| 
و 9و0 0 و0 باغ کوثر مین م0 078۳ و09 توبن زهرا ۳10 
سم سس وک 0 ی ی سرت سس ی دیص سس 9۱0112 ( 
۱0۳ 09 0.6 جیرکول 0 075 09 یت ۳10 
یی ال سس سس 08)88 12 9 
و 9و0 05 و0 0 0 051 088 ۰ ۳10 
کی تا ی ی وا میا سا 220۱ 0[ 
نیروگاه رجایی 
ب- ۱ 
و و0 075 090 تِِ ِ وین مهن 035 08 ِِ ۳10 
۱ ۳ 11۳0190 
وی و0 060 کوهین ا 0 072 و0 گزوين ۳10 
صتطاییک 022/18 


و ** به تربیت همبستگی معنی‌دار در سطوح ۵ و ۱ درصد را نشان می‌دهند. 
۲ 0610600 1 200 5 ۵۶ ولع7ع1 26 ممتاجاهمم امهه‌تنموزو لوح ۶ 


در جدول ۴ مشاهده می‌شود در همه ایستگاه‌ها به‌جز جیرکول 
میانگین دمای هوا بیشترین سهم و مهم‌ترین مولفه را در افزایش 
دقت تخمین 110 یا به‌عبارت‌دیگر بیشترین سهم را در کاهش خطای 
تخمین ۳10 دارد. ماجوزی و همکاران 
(2021 ,.۸1 ۵ ۷]۵[021) نتایج مشایهی را ارائه نمودند. هم‌چنین 
سرعت باد کمترین سهم را در افزایش دقت تخمین مقدار 10 دارد. 
در این تحقیق از رگرسیون بردار پشتیبان نوع دوم استفاده‌شده است. 
رگرسیون بردار پشتیبان دارای پارامترهای تنظیم ) و نا و 
پاامترهای هسته (۷) است. در این پژوهش نوع کرنل, از هسته تایع 
پایه شعاعی (۳) استفاده‌شده است که این هسته شامل پارامتر گاما 
() است. 

با توجه به جدول ۵ نتایج مقایسه سه روش در برآورد مقدار 
0 نشان داد که در هر دو مرحله آموزش و آزمون دقت روش 
۵5 نسبت به دو روش "1۳ و ٩۷۷‏ به‌طور نسبی بهتر بود. 

مقایسه میانگین نتایج حاصل از سه روش بر اساس آماره 
۱1/9 در ایستگاه‌های سینوپتیک و کلیماتولوژی استان قزوین در 
مرحله آزمون» در شکل ۲ ارائه‌شده است. در شکل ۲ مشاهده می‌شود 
در ایستگاه‌های سینوپتیک نتایج روش ۳ و 5۷۷1 مشابه و با نتایج 
روش ۸5 تفاوت معنی‌داری دارند. در ایستگاه‌های کلیم‌اتولوژی 
نتایج سه روش تفاوت معنی‌داری ندارند ولی میانگین خطای دو 


روش "18 و ۸((۷5 تقریباً مشابه و کمتر از روش 5۷ می‌باشند. 

نتایج نشان می‌دهد (جدول ۵) که میانگین خطای سه روش 
بهکاررفته در این پژوهش در بوئین‌زهرا قزوین و تاکستان که تقریبا 
ارتفاع یکسانی نسبت به سایر ایستگاه‌های سینویتیکی دارند. در هر دو 
مرحله آموزش و آزمون کمتر می‌باشد. افزایش دماء افزایش 1710 و 
تغییرات زمانی بارش در دشت‌ها باعث کاهش منابع آب و افزایش نیاز 
آبی گیاهان شده و در صورت تصمیم‌گیری نامناسب برای سازگاری با 
شرایط آتی باعث خسارات جبران‌ناپذیر خواهد شد ( ذ۷هنطک1 :۳۵2۵1 
0 ,.۸ 67). بنابراین با توجه به اینکه این سه منطقه به‌صورت 
دشت هستند و ميانگین دماء سرعت باد و 710 در این ایستگاه‌ها 
نسبت به ساير ایستگاه‌ها بالاتر می‌باشد (جدول ۱ بنابراین جهت 
تأمین نیاز آبی گیاهان این مناطق نیاز به تصمیم‌گیری مناسب جهت 
مدیریت منابع آب می‌باشد. در جدول ۱ مشاهده می‌شود» ضریب 
تغییرات دما برای ایستگاه‌های بوئین زهراه قزوین 9 تاکستان برابر ۶ 
و ضریب تغییرات 10 برای ایستگاه‌های بوئین‌زهراء قزوین و 
تاکستان به ترتیب ۰/۵۷ ۰/۶۰ و ۰/۶۱ می‌باشد. با توجه به یکسان 
بودن ضریب تغبیرات دما و 1770 و همچنین در ارتفاع تقریباً ی کسانی 
یکسان می‌باشد (شکل ۲). 


+ 29 تادز۵ن 204 نج وود ور مرگ صری صکوق یو مامتا نبا ربازک رات 
که جی ) بوخ یالیو . ) مج موم وب چم 6 62 ۳ ۳ آمو 


11 "0 18 ا 
۶1 ا 1 4 
1 0 ۱ ۱۱ 
۴ ا ۷ 6 
1۱ : 1 
+ 9 9 ۶۱ 
8 ا ۱۱ 0 
9 1" ۱1 4 
6 وا وا و 
۲ ‌# 6 ۲ 
۹ 0۱ ۱ ‌ 
زا ۲1 3 ۱ 
9 1 3 ا: 
۰ ۳ 6 ۱۱ 
11 3 ۷ 9 
۳۳ 0 اه 0۱۵ ۱0 و۱۹۱۱ ۵ | راخ را و۱۱ ۱۳ ۰ ۱۱20 
00۵ ۳۳ 0111 0 
پردسصو و در ورم و ار وود هر ۱۵۵۸ ۲03 ی تا 
1 ۲ ب هن ريت 6 و 62 


و "۱ 
,ونم و ۲7 روم اس رجعم ۳۲۶ رو 
٩3171 ۳ ۳35۳۲ ۳۵۱۱30۵ ٩۱۱۸۲۲۷ 0 ۲3۱۱۹۱۱313۷۵۳۷ -‏ ۳115 
پیب مت ای جع رود لین ۲6۳ کم کو + بر ۳6۲ 


۷۳۰ 


نشریه آب و خاک. جلد ۳۶. شماره ۶ بهمن - اسفند ۱۴۰۱ 


صداقت و همکاران. ارزیابی تبخیر تعرق مرجع با استفاده از روش‌های داده کاوی و مقایسه ... ۷۲۱ 


جدول ۴- اهمیت عامل به‌دست‌آمده با روش 1۳۳ 
۵ ۱۱۳ ها 0مصنه)ان 12600 6 0۶ ۵۲۲۵۳06۵ جرج -4 ما1 


میانگین رطوبت میانکین 


۳ ایستگاه‌ها 
آفتابی باد نسبی دما 
ص ۳۷۹ :۱ ی [ عصمناداه 
ایستگاه‌های کلیماتولوژی 
مصمتاقفاه زعمامتفصنلن 
البرز 
068 032 053 100 لبرز 
یاس 0018 ۳۹ 
095 034 02 100 کوتز 
سس سس سس 69988 طلعه ظ ی 
080 036 1.00 057 تن کول 
ی سس 60ص تست 
09 035 0:0 100 مکسلی 
سس سس سس 
نیروگاه رجایی 
099 043 00 1.00 نهز۳2 


با توجه به جدول ۵ و شکل ۲ دو ایستگاه رازمیان و سیردان با 
وه رام اد سطه فررای اه تیه سای ها 
سینوپتیک کمتر می‌باشد. اما در هر دو مرحله آموزش و آزمون 
3 آن‌ها بیشتر از ساير ایستگاه‌ها است. و درنتیجه دقت 
هواشناسی در دسترس در این دو ایستگاه (1-2192) محدود بوده 
است. از طرف دیگر با توجه به نتایج حاصل از تحقیقات شهریاری و 
همکاران (2019 ,.۵1 6 2۳772۲) که ایستگاه‌های سینوپتیک 
استان قزوین را از نظر عوامل اقلیمی بررسی کردند» مشاهده شد 
بیشترین درصد رطوبت نسبی (۶۰ درصد) را ایستگاه سیردان و 
کمترین درصد رطوبت نسبی را بوئین‌زهرا (۳۳ درصد) و تاکستان (۴۸ 
درصد) دارند و همچنین بیشترین سرعت متوسط باد سالانه را 
ایستگاه‌های سیردان و کوهین دارند. 

با توجه به جدول ۱ ضریب تغییرات سرعت باد در ایستگاه 
سیردان (۰/۶۵) بالا می‌باشد. بنابراین احتمالاً به علت شرایط اقلیمی 
ویژه‌ایی که در سیردان و رازمیان وجود دارد و همچنین محدود بودن 
می‌باشد (شکل ۳). مقادیر آماره‌های 13۳151 182 و :12۷181( روش 
۸۱۷9 در سیردان» به ترتیب برای هر دو مرحله آموزش و آزمون 
برابر با ۰۰/۴۳ ۰/۹۵ و ۰/۱۴ می‌باشد. درروش "1*۳ مقادیر آماره‌های 


ساعت سرعت میانگین رطوبت میانگین 


۳ ایستگاه‌ها 
آفتابی باد نسبی دما 
ص لک 6 ۶ ۱۰۸ زر عصمناداه 
ایستگاه‌های سینوپتیک 
عصمتاهاه متام0ظ5 
065 037 06۱ 100 آییک 
ی سس سس >8بوطاه 9 
030 034 040 100 اسماعیل آباد 
ید سس ۸084 086۵1 
آ 
057 019 058 1.00 ك 
یت یساس سس گت 
| 
053 018 02 100 بولین زهر 
ی سس سس 928]78تاظ "1 
043 057 031 1.00 ِ 
صعادی2] 
036 045 041 100 رازمیان 
ی یا مارا ابسسس لاس7 15211188 . 
۱ 
038 039 024 100 سیردان 
مخ و ۱ 
024 034 059 100 فزوین 
ی سس یت 
046 024 042 1.00 کوهین 
ی ی ی تا ی ی 
049 037 039 100 0 
طمردایصمله۷]۵( 


۴ 82 و ۱1۷/5 به ترتیب برای آموزش برابر با ۰/۴۵ 
۴ ۰/۱۵ و برای آزمون برابر با ۰۰/۴۷ ۰/۹۳ و ۰/۱۶ می‌باشد. 
درروش ٩۷‏ مقادیر آماره‌های 15 2 و ۱۷/5 به ترتیب 
برای آموزش برابر با ۰۰/۴۶ ۰/۹۴ و ۰/۱۵ و برای آزمون برابر با 
را در نظر بگیریم در هر دو مرحله آموزش و آزمون بیشترین دقت 
تخمین 10 در ایستگاه قزوین مشاهده می‌شود. در ایستگاه قزوین» 
درروش ۸۱5 مقادیر آماره‌های ۷/61 22 و ۱/6۳ به 
ترتیب برای آموزش برابر با ۰/۴۰ ۰/۹۷ و ۰۰/۱۱ و برای آزمون برابر 
با ۴۰« ۷« و ۰۱۱ می‌باشد. درروش 1۳ مقادیر آماره‌های 
۴ 82 و ۱1۷6۲ به ترتیب برای آموزش برابر با ۰/۴۲ 
۶ ۰/۱۲ و برای آزمون برابر با ۰۰/۴۲ ۰/۹۶ و ۰/۱۲ می‌باشد. 
درروش ٩۷‏ مقادیر آماره‌های 15 2 و ۱۷/5۳ به ترتیب 
برای آموزش برابر با ۰۰/۴۴ ۰/۹۶ و ۰/۱۲ و برای آزمون برابر با 
۴ ۰/۶ و ۰/۱۲ می‌باشد (جدول ۵). با توجه به شکل ۳ و ۴ 
مقادیر 118۷1610 در ایستگاه‌های کلیماتولوژی کمتر از سینوپتیک 
می‌باشته بنابراین دقت مین و۲۳ ایشگاه‌های کلیضاتولژی بالافر 


است. 


۷۲ 


نشریه آب و خاک. جلد ۰۳۶ شماره ۶ بهمن - اسفند ۱۴۰۱ 


جدول ۵- نتایج روش‌های مختلف در برآورد مقدار 171 در ایستگاه های هواشناسی 
کصماهاو آینع ۳6۵۲۵۵۲0۱0 عظ) صا ۴۳ 0۶ مصاجنای 0۲ موی ولمطامصه اصمه]ان0 01 دامع 1۳6 -5 12016 


آزمون (1650) آموزش (1۳210) 
ایستگاه‌های سینویتیک 
عصمتاهای 00 

(۷9۳ ۳2 1۷۳ ۱/۵ 0 1۳11۲ (۷9۳ 1 ۱/۱۰۵ ۱/۵ 1 1۳1۲ 

1 1۳ 036 0.96 ۰ .-02 012 099۵. _]6 039 095 0017 013  .9909۵2( ۰.۹4 

5 099۵ 012 0012 095 036 6( .099۵ 012 03- ۰ 0.96 036 5( 9 
0 ۱ 0 ۵9 و0 00 و09 040 ۰ 5۷/۸۲ 9[ 

۴.5 099۵ 013 90.-. ۰ 0.95 0352 6( 099۵ 013 5.-. ۰ 0.96 048 1۳ آری 

6 . 0099۵ 012 6.-. ۰ 0.96 047 6( .099۵ 012 5-. ۰ 0.96 047 ۸ ص 
 099‏ ق۵ 008 096 0:50 0 ۵9 و0 ۵00۵ و09 060 ۰ 5۷۸۲ اس و 

6 . 099۵ 012 0-. ۰ 0.96 044 6 : 099۵ 012 ۰-6 ۰ 0.96 042 1۳ اتاخیل آراد 

6 0099۵ 011 0- ۰ 0.97 041 6( .099۵ 012 2.-. ۰ 0.96 042 5( 9 
099 ق006۵ 096 047 ار یی 60 دم 

9 1۳ 044 0.96 ۰ ۰-6 011 099۵. (6 044 0.96 ۰ .-05 0.12 0099۵ ۳.6 

6 . 099۵ 011 0.000 096 045 6]: 099۵ 011 0- ۰ 0.96 045 5( ات زهر 
099 ۵2 00 096 047 ا کف0 ۵9 مه0 ف008 06 08 کر 

۳.6 099۵ 013 0001 096 048 6]_ .099۵ 012 03- ۰ 0.96 044 1۳ تاکسا 

6 099۵ 011 0003 0917 043 :]099۵ 011 0000 096 043 5( تب 
6 09 ۵2 004 096 046 ا ۵۵6 ۵99 ۵ 00۵8 0۵6 ک0 ۰ 5۷/۸6 ۱ 

1 1۳ 045 0.95 ۰ .-0 013 099۵ (5 046 0.95 ۰ .-8 0 14 099۵ : 5 

7 ۸ 047 0.95 ۰ .-2 0 14 099۵ _(5 047 0.95 ۰ .-45 0 14 099۵ (0.5 


4 098 014 8 095 0.46 5 099 014 0005 095 043 ۳۳ معلم کلایه 
95 0.99 ۰.3 0024 095 042 95 099 ...013 0008 کو0 ...042 ...متخ( 
5 099 ۰ 014 8 095 044 4 098 . 014 008۵ ک 09‏ کف0 و طهوهایملد۷۵ 
ایستگاه‌های کلیما تولوژی 
عممناهاه رعماماهصنن 
۳۳ ۳۵8 1۳28 5۳ ۰ 32۳6۲ ۳۳ ۵۵8 1۳88 ۳ ۰ 32۳8۲ 
6 099 013 0012 096 042 1 097 099 012 0002 097 039 ۱.۵ مار 
6 099 3 0014 096 0.42 6 099 .۰ 012 . 0015 096 04 ولال۸ 
06 0.99 ار 04 0043 09 046 6 0.99 013 ۰.35 0.95 ۰ 0.45 9۷ آععع۱ 
7 099 ۰ 0.10 1 097 0.35 7 099 009 ۰ 0001 097 033 .10 ۲ 
7 099 ۰ 0.10 ۸ 097 .۰ 0.35 7 099 .۰ 0.10 0 097 035 ول((۸ مت 
07 0.99 009 ۵9 038 ۵96 ۵.99 ال 0 فا ی ۱ ۳ ۲۵۲۵ ۳۵۳۶۱ 
۰6 099 012 2 096 044 6 099 011 0.1 7 0۵ دا ار 1 نِِ 
7 099 011 0 097 .۰ 0.40 7و 9و0 1 0009 097 .۰ 039 _._ ول((۸ ف 
6 0.99 و 004 ۵9 046 م۵ و ۵۵9 م09 م04 5۷۲ ی 
۰6 099 013 3 096 044 7 099 011 03 097 ۰ 038 ۹3 با کوئر ‏ 
6 099 013 0 096 . 0.44 6 099 .۰ 013 0014 096 .۰ 042 ولال۸ جع ور 
 8 ۵48 ۵۵ 2‏ م09 ق04 5۷۲ ی 
5 099 ۰ 014 .۰ 0012 095 044 6 099 013 .۰ 0002 096 043 . 3 چیرکول نك" 
95 099 ۰ 0.14 4 095 0.46 5 099 . 014 .۰ 0007 095 046 ول( امن 
4: 0:9۵ 015 4 094 049 4 099 . 015 . 007 094 049 ۰ 5۷۷۲ 


صداقت و همکاران ارزیابی تبخیر تعرق مرجع با استفاده از روش‌های داده‌کاوی و مقایسه ... ۷۳۳ 
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شکل ۲- مقایسه :۱1*151 سه روش در برآورد 1710 در مرحله آزمون 
5عاز اجع) ۲۳6 صا و۳ ۵۶ مدتاههایی 1۵۲ مل‌مطاعصهر معط ان ۱۱۱۳۱۲۹۲ م6 «معزنهم‌همن -2 ۲12۱۳۲6 


نتایج نشان می‌دهد (جدول ۵) که دقت سه روش به‌کاررفته در 
آکت تقوفتن فن ایسگاه فیروگاه رای شب سای آسمااضاین 
کلیماتولوژی, در هر دو مرحله آموزش و آزمون بیشتر می‌باشد و 
۱8 روش‌های بکار رفته در این ایستگاه نسبت به بقیه کمتر 
می‌باشد. و همچنین دقت تخمین در ایستگاه جیرکول نسبت به سایر 
ایستگاه‌های کلیماتولوژی در هر دو مرحله آموزش و آزمون کمتر و 
میانگین خطای بیشتری دارد. افزایش ارتفاع سبب کاهش دما و به‌تبع 
آن کاهش مقدار 710 و افزايش رطوبت نسبی می‌شود ( 0 ۵216 
5 ,./0). در ایستگاه جیرکول به علت افزایش ارتفاع و کاهش 
مقدار 0 کاهش دقت تخمین در اين ایستگاه منطقی می‌باشد. اما 
در ایستگاه تیروگاه رای کهبا افزایش ازضاع دقتت مین 
افزایش‌یافته است نیاز به بررسی بیشتری می‌باشد. در این منطقه به 
علت فعالیت نیروگاه رجایی و تولید فراوان آلاینده‌های مضر که از 
ترسیب آن‌ها موجب آلودگی آب سطحی, خاک و گیاه می‌شود. حضور 
این آلاینده‌ها باعث باران اسیدی و درنتیجه نابودی پوشش‌های 
گیاهی و ساختار خاک می‌شود (2011 .۸ 21 ن2). در این منطقه 
باغات و مزارع کشاورزی در جنوب» جنوب شرقی» جنوب غربی این 
نیروگاه قرار دارنده از طرفی جهت باد غالب در استان قزوین به سمت 
جنوب شرقی بوده (2019 ,./۵ 61 9۳۵72۲) و بیشترین میزان 
ترسیب آلاینده‌ها در این جهت می‌باشد که باعث افزایش ۳۲۲0 و 
کاهش رطوبت نسبی و همچنین خسارت به پوشش گیاهی می‌شود. 
این مسئله در جدول ۱ مشاهده می‌شود, میانگین رطوبت نسبی در 
آشتگاه تبرزگاه رجایی.تا مجند ازفاع بالاش کمعر از سایر ایستگانها 
است. کیخسروی و همکاران (2014 ,.۵ ۵ 161100972۷1) در 
اییشگاه‌های ویک اسان کرماشاه واظ را مخاسية کرونته و 
نشان دادند که عامل ارتفاع باعث اختلاف در عناصر اقلیمی و به‌تبع 
آن مقادیر 1710 شده است. 

در پژوهش حاضر با افزایش ارتفاع از ۱۳۰۰ متر به بالاء میانگین 


خطا به‌طور قابل‌توجهی افزایش می‌یابد (شکل ۳). علت آن‌هم احتمالا 
به دلیل تغییرات دمای هوا با افزايش رطوبت و ۳ آن بر مقدار ۳۴0 
افزايش عرض جغرافیایی از میزان تشعشع ورودی به زمین کاسته 
می‌شود. بنابراین انتظار می‌رود که مقدار 1210 در عرض‌های جغرافیایی 
بالات کمتر از عرض‌های جغرافیایی پایین‌تر به دست آید ازجمله 
دلایل اين امر می‌تواند مرتبط با در نظر گرفتن عواملی مانند عرض 
سایر متغیرهای آب و هوایی دست‌خوش تغییرات قرار می‌گیرد (۷۵2 
7 ,.۵ 6). بنابراین نتایج تحقیق حاضر نشان می‌دهند که اثر 
نامشخص است. کنگ و همکاران (2021 ,.4 6 ع20>) نتایج 
مشابه‌ایی ارائه کردند. 

با توجه به جدول ۵ مقادیر آماره ۷1۳817 در مرحله آموزش و 
آزمون اکتا مقادیر کوچک و نزدیک به صفر می‌باشد بنابراین با توجه 
به کم بودن مقادیر آماره 2۷187 نمی‌توان گفت مدل کم‌برآوردی یا 
بیش برآوردی قابل‌توجهی داشته است و درنتیجه می‌توان از آن 
چشم‌پوشی کرد (2020 ,.21 6 81721). در جدول۵ مشاهده می‌شود 
۴ می‌باشد و هم‌چنین شاخص 4 هم در همه ایستگاه‌ها در هر دو 
می‌باشد. از طرفی فاکتورهای مختلفی بر عملکرد روش‌های 
یادگیری ماشین مانند الگوریتم آموزش ۸۱5 یا تنظیمات فرا 
پارامترهای روش‌های مورداستفاده. دخیل هستند ( ,.۵1 6 ۳۲۵۲۲۵۲۵ 
9 بنابراین در این تحقیق انتخاب عوامل موّثر در تخمین مقدار 
110 مناسب بوده است که مقادیر آماره‌های ۳۳ و9 و بالا بوده آنبتت: 
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شکل 6- تغیبرات خطا با افزا 
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نتیجه گیری هم‌چنین به‌شدت تحت تأثیر شرایط منطقه‌ای مانند توپوگرافی و 
ارتفاع می‌باشد. نتایج این تحقیق علاوه بر بررسی بهتر 1010 به 
2 شناخت عوامل تأثیرگذارتر در هر منطقه کمک فراوان می‌کند و 
ی ی ام رسای بازهي ها ماه ارآن ناه کردمو وتا 
دارد. در این پژوهش دقت تخمین مقدار 1 با دو روش 
۸5 و 6۳ نزدیک به هم و نسبت به روش 5۷1۷ بالاتر بوده 
است. نتایج نشان داد مقدار ۳0 نه‌تنها بر اساس دمای هوا بلکه 


ممکن است تحت تأثیر عامل‌های دیگر مانند آلودگی هوا تغییر کند» 


نتایج پژوهش حاضر نشان داد مقادیر ۴10 سامانه نیاز آب» 


در منطقه مطالعاتی حاضر مشخص شد که میانگین دمای هوا بیشتر 
از سایر پارمترهای اقلیمی بر 10 تأثیر می‌گذارد. با درک این مستئله 
چنانجه دمای هوای یک منطقه زیاد شود باعث افزایش شدید در 


۶ شریه آب و خاک. جلد ۳۶. شماره ۶ بهمن - اسفند ۱۴۰۱ 


مقدار 70 و متعاقباً افزایش نیاز آبی گیاهان خواهد شد. توجه به این نقش ارتفاع و عرض جغرافیایی است. نتیحه نهایی اينکه الگوریتم‌های 
نکته خواهد توانست از بروز تنش آبی و کاهش احتمالی محصول مختلف روش‌های داده‌ کاوی» در شرایطی منطقه‌ایی مانند سیردان و 
جلوگیری کند. با توجه به متفاوت بودن ارتفاع و عرض جنرافیایی رازمیان که دادههای اقلیمی محدودی در دسترس می‌باشده می‌توانند 
ایستگاه‌هاء و از طرفی شرایط اقلیمی متفاوت در این ایستگاه‌هاء در مطالعات آینده برای توسعه مدل‌های تخمینی 1710 مناسب باشند. 


می‌توان چنین نتیجه گرفت که نقش اقلیم در تخمین ۳۲0 مهم‌تر از 
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05://001.0۲8/10.1007/5001704-019-03007-۰ .1459-1477 :139 چوماها۱0۱۵) :م۸ 


6۱۱0و 20 2710 طا ولو لقصه 0صقط امنهر (داصصم .(ر2005) .ظ رصهکفطاوهو عک ررظ مطمافجته تفع لطومصونا ریز رلملمفگ .. 


۳۵۲۹12(۰ 12) 73-81 :54 م۷۷۵۵ ۵۵ ۲۷۷۵۸۵۲ .صهع۱ معاعده 20 لهاجهع ۵1 وهم1عع۲ 


مط۱ طهزصدط عطا که مصتاه0مص مه ممتامتایی ممتامام‌مهامم۵۷۵ ۴۵16۲۵888 ,(201۱9) .9 رصلعل ی بنکا.۷ رلعههه . 


38-۰ :156 ...0۵و ...0 عع80100ع .0۱94 ...۵۴۵96۲ بعصتصنجع1 ...4660 هصتویا..۰ 18012 
0۰ 6/1 10.101 /1005://001.0۲2 

۵ 0027۷1 طا قط۵97) فص فطع گم مصتاهع1۵ لقصصتامن) 201۱9(۰) یک رتاهممصصعطوها۲ ی ریخ رکقلطمصنهت بیبط رطف 
.(619120ظ . ظ) ,19-34 :29 . ممهاظ. امه۵۵ا کصه واموهنوهعع .. .ععامصصحتهم .. متاقصتاه .. صم .۰ 60ق2ظ 
0۰ 08/260 10.221 /2 001.۰0۲ //:۱)۲05۹ 

5 0 111۷ هه عصاممرصظ ز10 وملوماه0ه مصتااناوه ما۲2 ,(201۱9) .0 روهتم0 ی ر.و باتک روط متاله۷ رل بلط 
1 48۳۱۵۱۸۵۱۱۱۲۵ ۱ 200 ۵۱۱ و۵۵۵ .عجمتاهاه تقانصته عجمصصح ممتامصصتاوه طملاهتامقصه ۵۷3۵0۲۲ ۲616۲۵8۵۵ ۲0۲ 
219.۰ 10۰1016/10 /0005://001.0182 .70-81 

۲ کلام تامتاومن . تقوم عط1۵۷۵10۵ (02020 ۲۷۰ رتعصتناظ عگ 1.۹ ما62 ,۵ ۳۱22۷6 بظ بمطع19 
٩011 26۲‏ 1۳6 ,۲ :ظ۷۹)۵و ۳۲۴۵۳ ۴۱۷ عط 4صه 0مط)عصه ممتاهتممهبه مط طا 0عسعدعصه وللمو لمصمتاقصعاص تم فحمتاع‌صیظ 
.۰ 2 .3425۰ :12 ۲/۵۸۵۲ 09۲۷6۰ عمتاصهاع۲ 

الب مج 0عاتصنا رلعکصته: 10۲ املمهصه 2۵تهص-ترهاوظ ۵ اصمصووعووه ممصممص۲۵۲]0 2020(۰) .9 ملق عک ریک ما20 
16/[۰2820726.2020۰ 10.10 /۸۵5://01:0۲82 106402۰ :242۰ ۲/۵۱۱۵ ۲۷۷۵۸۵۴ ۸۵۲۶۵۲۵ .عصمتاتصمی ۵26109 1۳11 

مه ۷۵۵۵۲ 0۲ممتاو ۵ رالانماهء‌نآنود ,(2013) ۲۰ رممهاه1 ۲۱۵95۵182۵0 6 ره ملطه۲2 و ,2عظ/۷2 وبا بلتهطاه 1 
1 562168 ۱۵۵۱۵۷ صمتاهنکصدام هه وم مقامم عصتام0مصظ عم مدرد ممصععع 22نظمنافه 202001۷6 220 
.1-012-0332-۰ 10۰1001750027 /1005://001.012 :575-588 

عنلفهنل وه 0ععهه همناتلعتن صمتامنم‌نصماممد۱ه معمعتعاهد وانظ۳ 2020(۰) ۱۷ بلقصما۱ ک بم۷.2 بهیزفک و۲ موعصتاً 
0 01 160۳61۱6۵ .وطماهناوه لدمنترجصه مهد فمیوتصطعا ‏ مصتصنص حاع . اصمتمله ‏ راومه فعمت‌ههع 
605://001۰01۲82/10۰10017/500704-020-03225-۰ ۰763-773۰ :141 روماهمامصنت 

0 04 6181۱6۵۲۱۵ ۵۲ 05اه 6 10 ۱60۳۱۱۵ آ6عا5 51۵ 0 ۵۳۸۵۲۵ 11:2 .(2000) ۷۰۲۰ کللط20 ۷ 
.21:1003-۰ ,۷0۲16 ۱۱6۷ ,50۳16۲ ,56161166 

لمطمتمصعهاما ,ععصممموم املمصه 0 »م10 0مصلله۲ ۸ ,(2012) کل هتفه۷ عک .9 مبط0ععطامط ر.ژ.) ملامطالز ۷۷ 
۰ ۵۰۲۵/10 .2088-2094 :32 7«ع2)010ص011 ۵۶ اهرتتاه [ 

۵ 1۳002015 باتااماه۷2 متفصتاه عصتاه‌عاه1۳ 2017(۰) .2 مط ۷‏ وک م۷۷ ول یلا ون) بنا .60 ب۲۵98 وبا ۷218 
596-۰ :48 عم سوماه۳۵۲ فصن ۱۷۷ رصتفدظ بع۷ ۳ مطمه1. فطا طا. صمتاهتاممممتممده . تقتطامع . 2۳00 
۰ 0 0۹:0( 
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